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Abstract 

Dirac hamiltonian on the Poincare disk in the presence of an Aharonov-Bohm flux and a uniform 
magnetic field admits a one-parameter family of self-adjoint extensions. We determine the spectrum and 
calculate the resolvent for each element of this family. Explicit expressions for Green functions are then 
used to find Fredholm determinant representations for the tau function of the Dirac operator with two 
branch points on the Poincare disk. Isomonodromic deformation theory for the Dirac equation relates 
this tau function to a one-parameter class of solutions of the Painleve VI equation with 7 = 0. We analyze 
long distance behaviour of the tau function, as well as the asymptotics of the corresponding Painleve VI 
transcendents as s ^ 1. Considering the limit of flat space, we also obtain a class of solutions of the 
Painleve V equation with (3 — 0. 
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1 Introduction 



It has been known since [111 |34] that the two-point correlation function of the 2D Ising model in the scaling 
limit is expressible in terms of a solution of a Painleve III equation. This remarkable result turned out to 
be a special case of a more general phenomena, described by Sato-Miwa-Jimbo (SMJ) theory of holonomic 
quantum fields and monodromy preserving deformations of the Dirac equation One of the central 

objects in SMJ theory is the r-function of the Dirac operator acting on a suitable class of multivalued 
functions on the Euclidean plane. 

The SMJ r-function admits a geometric interpretation [27l[28]. Loosely speaking, it can be obtained by 
trivializing the det*-bundle over an infinite-dimensional grassmannian, composed of the spaces of boundary 
values of certain local solutions of the Dirac equation, where different points of the grassmannian correspond 
to different positions of the branch points on the plane. The same idea was earlier used in [26] to show that 
the T-function of the Schlesinger system can be interpreted as a determinant of a singular Cauchy-Riemann 
operator. A simple finite-dimensional example of this construction arises in the study of a one-dimensional 
Laplacian with (5-interactions 20J. The most important thing about the geometric picture is that it allows 
to find an explicit representation of the SMJ r-function in terms of a Fredholm determinant, thereby giving 
a solution of the deformation equations. 

The simplest way to generalize the above setup is to replace the plane by an infinite cylinder. In this case, 
the deformation equations and a Fredholm determinant representation for the r-function of the corresponding 
Dirac operator were obtained in [19]. These results provide a shortcut derivation of the PDEs satisfied by 
the scaled Ising correlation functions on the cylinder [18 and of the exact expressions for the one- and 
two-particle finite- volume form factors of the Ising spin and disorder field [H jSj [13] (and, more generally, of 
twist fields [To]). 

In 29J, Palmer, Beatty and Tracy (PBT) have extended the SMJ analysis of isomonodromic deformations 
to the case of a Dirac operator on the Poincarc disk (see also earlier works [25], [30] on this subject). The 
associated r-function in the simplest nontrivial case of two branch points was shown to be related [291 to a 
solution of the Painleve VI equation 



cPw 1/1 1 1 \ f dwV [\ 1 \ \ dw 



ds^ 2 \w w — \ w — s J \ ds J \s s— 1 w — s J ds 

w{w — 1){'W ^ s) f ^ [3s ^ s — 1 s{s — 1) 



(1.1) 



5^(5—1)2 \ {w ~ 1)'^ {w — s)^ 

with only one fixed parameter. Before stating the PBT result in more detail, it is useful to reformulate it 
in a slightly different way. The Dirac operator, considered in [29], is simply related to the hamiltonian of 
a massive Dirac particle moving on the Poincare disk in the superposition of a uniform magnetic field B 
and the field of two non-integer Aharonov-Bohm (AB) fiuxes $1,2 — 27r^i^2 located at the points ai and 02. 
Without any loss of generality, one may choose — 1 < i^i^2 < 0. It is preferable to work with the hamiltonian, 
as it is a symmetric operator that can be made self-adjoint after a proper specification of the domain, and 
many assertions of [29] (e. g. symmetry of the Green function, non-existence of certain global solutions of 
the Dirac equation, etc.) immediately follow from the self-adjointness. 

Write the disk curvature as — 4/i?^, denote by m and E the particle mass and energy, and introduce two 

dimensionless parameters b = /i = Vt'^' +'^^ ^ ^ turns out that the r-function associated to 

the above hamiltonian depends only on the geodesic distance c?(ai, 02) between the points ai and 02. If we 
further introduce s = tanh^ d{aua2) ^ ^.j^gj^ ^^^j-^ |-,g expressed [15] in terms of a solution of the Painleve VI 
equation (jl.ip 

— lnr(s)= — + (1.2) 

where X ~ 1/2 — 1^1, A = 2 + z^i + 1^2 — 26, and the values of the PVI parameters are given by 

A2 (A -1)2 1-4^^2 
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Actually, the paper [29j is concerned with the case E — (the mass term in the Dirac operator is not of the 
most general form) . We include this parameter from the very beginning because it will be shown below that 
the final answer for the r-function depends on E only via the variable fj,. 

The aim of the present study is to solve the remaining part of the problem — that is, to compute the 
PBT r-function and to investigate its asymptotic behaviour, which can be used to specify the appropriate 
initial conditions for the equation We summarize our results in the following theorem: 

Theorem 1.1 The PBT tau function admits Fredholm determinant representation 

T(s)-det(l-L,,,,L',^,), (1.4) 
where the kernels of integral operators L^^g, L'^ ^ are 

K,siP,(l) = Ly.s{-p,-q), 
Is = arctanhy^ = d{aya2) g g M. The functions p{p), J^i^ip, q) are given by 

22Mp(i_,_2^) 



pip) 



r (/i + i + f ) r (/i ^ i - ^ 



2 ' 2 / V' ' 2 2 , 
tt/2 7r/2 



^uiP,q) = n" i de I dx I dy 



27r2 J J J ^ (e^^^n^-y) + l)(2 + 2cosht 

— oo \ / ^ 

X (sin x) ^+ ~ 5 (cos x)^~~f~^ (sin y) ^ ~ ^ (cos y ) ^+ ^ ~ 5 . 

To leading order the long-distance (s ^ 1) asymptotics of t{s) is 

1 - t(s) ~ Ar{l - s)i+2^ + O ((1 - sf+^^') ass^l, 

where the coefficient A^. is given by 

sinTTi/i sin7ri^2 r(/i -f 2 + — b)T(fi — vi + b)T(^ + 2 + — b)T(fi — V2 + b) ,^ ^, 

At = T, 9 . (1.5) 

[r(2 + 2^)]' 

A few remarks are in order. The formula (jl.Sp implies that the asymptotic behaviour of the corresponding 
PVI transcendent for /x > 1/2 is 

1 - w{s) ~ A[l - s)i+2^ + O ((1 - s)2+2a') as s ^ 1, 

A - (1 + 2^)^ A 

(Ai+ 1 + 1^1 - 6)(/i + 1 + z/2 - 6) ^ 

The transformation i/i i— s- i/i + c, 1^2 > ^^2 + c, 6 i-^ + c, /i i— > /i does not change the values (|1.3p of the PVI 
parameters a, f3, 7, S. However, our solution nontrivially depends on all four variables vi, V2, b, fi as can 
be seen from the above asymptotics. Thus we have constructed a one-parameter family of solutions of the 
Painleve VI equation with 7 = 0. 

We also note that in the case 6 = the PBT r-function was conjectured to coincide with the correlation 
function (01+1,^(01)01+1^2(02)) of U{1) twist fields in the theory of free massive Dirac fermions on the 
Poincare disk [J (particular case of Ising monodromy was later studied in more detail in [5] l^)- The 
asymptotics of r(s) as s ^ then follows from the known flat-space OPEs for twist fields. Long-distance 
behaviour is determined by a form factor expansion of (Oi+,yj (01)01+,^^ (02)). The method of angular 
quantization, employed in [7] for the calculation of form factors, does not seem to work quite well, as it leads 
to formally divergent expressions. A sensible answer for the infrared asymptotics was nevertheless extracted 
from them after a number of regularization procedures. The formula (jl.Sp . specialized to the case 6 = 0, 
proves the latter result. 
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The main technical problem arising in the direct computation of t(s) is the unknown formula for the 
Green function of the Dirac hamiltonian on the disk in the presence of a uniform magnetic field and one 
AB vortex. Such hamiltonian can always be made commuting with the angular momentum operator by a 
suitable choice of the gauge. Partial Green functions are then calculated relatively easily in each channel with 
fixed angular momentum; the difficult part is the summation of these partial contributions to a closed-form 
expression. We have solved this problem by writing radial solutions of the Dirac equation as Sommerfeld- 
type superpositions of horocyclic waves, similarly to a simpler scalar case [21] . This allows to perform the 
summation and to obtain a simple integral representation for the one- vortex Green function. 

This paper is organized as follows. In Section 2, after introducing basic notation, we describe the solutions 
of the radial Dirac equation and compute radial Green functions. Spectrum, self-adjointness and admissible 
boundary conditions for the full Dirac hamiltonian are also briefly analyzed. Finally, we write contour 
integral representations for the radial solutions and obtain a compact formula for the one- vortex resolvent. 
The PBT T-function is studied in Section 3. We start by giving a general definition of the r-function in 
terms of the projections on some boundary spaces. Next we introduce coordinates in these spaces, using 
the solutions of the Dirac equation on the Poincare strip. Explicit formulas for the projections in these 
coordinates are obtained in Subsection 13 . 31 bv analyzing the asymptotics of the Green function, computed in 
Section 2. These formulas give the kernels of integral operators in the Fredholm determinant representation of 
r(s). In Subsection [331 we recall the relation of the PBT r-function to Painleve VI equation. Subsection 13. 51 
deals with the derivation of the long-distance asymptotics of t(s). The analogs of the above results in the 
limit of flat space, where Painleve VI equation transforms into Painleve V, are established in Section 4. The 
Appendix contains a proof of the fact that the r-function depends only on the geodesic distance. 



^^%W = ^r^' ](^SU{1,1). (2.2) 



2 One-vortex Dirac hamiltonian on the Poincare disk 
2.1 Preliminaries 

Let us first establish our notations. We denote by D the unit disk |zp < 1 in the complex z-plane, endowed 
with the Poincare metric 

els'" =g,,dzdz^R^-^^^^. (2.1) 

(1 - |z|2) 

This metric has a constant negative Gaussian curvature — 4/i?^ and is invariant with respect to the natural 
SU{1, l)-action on D: 

az + (3 _ f a (3 

Pz + a' ^~ \ P a 

The hamiltonian of a Dirac particle of unit charge moving on the Poincare disk in an external magnetic 
field has the form 

where the operator K and its formal adjoint K* are given by 

K = \2D, + l-dAng^A , (2.4) 



K* = ^ \2D, + i d, \ng,, } , (2.5) 



9zz { ' 2 

and = dz + iA^, Dg ~ dg + iAg denote the covariant derivatives. 

Connection 1-form A — dz + Ag dz, which is considered in the present section, consists of two parts. 
Namely, we put A A^^^ + A''''\ where 

AiB) _ JBR^zdz-^zdz 

- 4 l-|z|2 ' ^^■^> 

= -'Ji('^-1^^4]. (2.7) 
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The first contribution describes a uniform magnetic field of intensity i?, since dA^^^ is proportional to the 
volume form dii = — g^z dzAdz. The second part corresponds to the vector potential of an AB flux $ = 27rz/, 

situated at the disk center. 

Introducing polar coordinates z = re'^ , z = re~^^ , one can explicitly rewrite the operators K and K* as 
follows: 



K = 
K* = 



R 



(1- 






r 




(1- 






r 





(l + 26)r 



(2.8) 
(2.9) 



BR 

Here, we have introduced a dimensionless parameter b — — — 
the disk curvature. 

2.2 Radial hamiltonians and self-adjointness 

Since the formal hamiltonian 



characterizing the ratio of magnetic field and 



angular momentum operator L 



corresponding to the vector potential + A^"'' , commutes with the 



— cfz , we will attempt to diagonalize them simultaneously. The 

eigenvalues of L are half-integer numbers + 2 (^0 G ^) and the appropriate eigenspaces are spanned by 
the spinors of the form 

, , / i«i„,i(r)e»''"^ 



The action of H leaves these eigenspaces invariant. One has 



u'/o,i(r) 
w/„,2(r) 



lo+v 



w;o,i(r) 



Hi = R- 



mR Ki 



Kr 



-mR 



where the operators Ki and are explicitly given by 



Ki = (l-r2)(9, + i±l) +(H-2&)r, 



r 

= -(l-r2)(a.-i-)-(l-26). 



Let us make several remarks concerning the solutions of the radial Dirac equation 

(Hi - E)wiir) = . 



(2.10) 

(2.11) 
(2.12) 

(2.13) 



It will be assumed that E G C\ ^{—00, — m] U [to, 00)^ and / is an arbitrary real parameter. We also introduce 
for further convenience the following quantities: 



S2)i?2 4, 4^2 



m + Ej VM - ^ 



All fractional powers in these formulas are defined so that fi and C± are real and positive for real values of 
E satisfying \E\ < to. 

Let us first look at the space of solutions of (|2.13p on the open unit interval / = (0, 1), which are square 



integrable in the vicinity of the point r = 1 with respect to the measure djjr = R 



rdr 
(l-r2)2 



, induced by the 
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Poincare metric. It is a simple matter to check that for any I G K. this space is one-dimensional (i. e. the 
singular point r = 1 is of the limit point type) and is generated by the function 



Wi (r) = 



^JJ^ r{fi - 6)r(M + 6) 



2/i 

X (l-r^)- 



Ti2^i) 



-1 ^-i 



2F1 - 6 + 1, ^ + & - /, 1 + 2^, 1 - 
:Fi (/i - &, /i + 6 - 1 + 2/i, 1 - 



(2.14) 



y^?^ V{pi - 6)r(M + 6) 



2/x 

X (l-r^)- 



r(2/.) 



-1 J 



2F1 (/i + 6, /i - 6 + 1 + /, 1 + 2/i, 1 - r2 



C+r'+i si^i (/^ + + 1, M - ^ + 1 + ^ 1 + 2/i, 1 - 



(2.15) 



If I e (—00, —1] U [0, 00), the limit point case is also realized at r = 0. The solution that satisfies the condition 
of square integrability in the vicinity of r = can be written as 



I 



(r) 



f T{^I-b + l + l) 
r(i + or(/i-6 + i) 



V 







r(/i-b + i + o 
r(2 + or(At-6) / 



\ C+r'+i 2Fi{^i + b+l,^i-b+l + l,2 + l,r^) 



X (1-r^) 

/ Ti^i + b-l) 



(2.16) 



(11 -) 



(r) 



r(i -or(/i + 6) 






r{pL + b-i) 
'^{-l)^i^I + b + l) ) 



X (i-r^)- 



C4 



2F1 {[L-b,\L^b-l, ~l, r^) 



(2.17) 



where the first formula corresponds to the case I > and the second one to I < —1. For / e (—1,0) both 
solutions (|2.16p and (|2.17p are square integrable at r = (the limit circle case). 

The functions Wi^\r) and Wi^^'~^\r) are linearly independent for I > —1, and the functions wl'^\r) and 
\r) are linearly independent for I < 0. One can show this, for example, by computing the determinant 



of the fundamental matrix constructed from these solutions: 

det (z^(^)(.),z^("-±)(r)) 



1 - 



(2.18) 



This implies that for E e C\ (^(— 00, — rn] U [to, cxj) j and I G (— 00,— 1] U [0, cxd) the equation (|2.13p has no 

square integrable solutions on the whole interval /. For I G (—1,0), however, there is a one-dimensional 
space of such solutions, generated by the function (|2.14p - (|2.15p . 

We now examine the issue of self-adjointness of the operators Hi . The above remarks can be summarized 
as follows: 

Proposition 2.1 Let us restrict the domain of the formal radial hamiltonian Hi, defined by Ii2.10\) - l2.1'2\) . 
to smooth functions with compact support in I . Then 

• Hi is essentially self-adjoint for I G (—00, —1] U [0, 00); 

• for I G (—1,0) the operator Hi has deficiency indices (1,1) and admits a one-parameter family of 
self-adjoint extensions (SAEs). 
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Assume that I G (—1,0). Deficiency subspaces IC± — ker \^Hj =F imj are generated by the elements 

w±{r) = wl^'^r) . (2.19) 

E—±im 

Different SAEs H^"^^ are in one-to-one correspondence with the isometries between /C+ and /C_. They may 
be labeled by a parameter 7 g [0, 27r) and characterized by the domains 



domi?/'^^ = |wo + c (w+ + e*^u;_) jwo G domi?,, c G c| . (2.20) 



It is also conventional to characterize the functions from the domain of the closure of Hi^'' by their asymptotic 
behaviour near the point r = 0. Namely, it follows from ((TTS)) and ^^J^ - ^^ that for w G dom (^H^''^ 
one should have 

lim cos I — H I (mRr) '' Wi{r) = — lim sin [ — H ) {mRr)^'^'' W2{r). (2.21) 

r^o y 2 4 y r^o \2 4 / 

Here, we have introduced instead of 7 a new self-adjoint extension parameter O G [0, 27r), defined by 

'e 7r\ _ 2-' r(-o 

*'"'2+4j - tan(I-i)-lf(TTO^ 



r(/i + b + i)r(M-6 + ^ + i) / A-^ /i^A ' " 
r(A-5 + i)r(M + fe-0 V/'^ + ^VV2y ' ^ 

where /x = = ^ ■ JNote that the choice B = — (— — ) is equivalent to requiring the 

regularity of the lower (resp. upper) component of the Dirac spinor at r = 0. 

Let us now consider full Dirac hamiltonian. Since the shift of the AB flux by any integer number is 
equivalent to a unitary transformation of H, hereafter we will assume that —1 < ly < 0. 

Proposition 2.2 Suppose that domH ^ C^f (D\{0}). Then 

• for i> = the operator H is essentially self-adjoint; 

• forv G (—1,0) it has deficiency indices (1,1) and admits a one-parameter family of SAEs, henceforth 
denoted by H^"'\ which correspond to those of the radial mode with Iq = 0. 

2.3 Radial Green functions 

Let us begin with the case I G (— 00,— 1] U [0,oo). Green function G£;.;(r, r') of the radial hamiltonian Hi 
can be viewed as the solution of the equation 

(^Hi{r) - e) GEAr, r') = 5{r - /) I2, (2.23) 

which is square integrable in the vicinity of the boundary points r = and r = 1. Standard ansatz 

At, w\"'^\r) ® ( w\^\r')Y for < r < r' < 1, 
GE.M-i I, , \' {t (2.24) 



4± 



I w^/\r) (g) for < r' < r < 1 



solves ()2.23p for r ^ r' and meets the requirements of square integrability. The sign "-f" ("— ") in the above 
formula should be chosen for / > (resp. / < —1). Prescribed singular behaviour of the Green function at 
the point r = r' is equivalent to the condition 

Gi5,;(r-f-0,r)-G£;,,(r-0,r) = -^^-y— ■ 
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Substituting p.24p into the last relation, we may rewrite it as follows 

1 _ J.2 



A^ydeiw(w\'\rlw\"^^\r) 
Finally, using (|2.18[) one finds that 



Rr 



Now suppose that I £ {—1,0). In order to find the resolvent of the radial hamiltonian hI^\ we need a 

solution of the radial Dirac equation (-ff/'''' — E)w'f '' — 0, which satisfies the boundary condition (|2.2f p at 
the point r = (square integrability near the point r = f is not required). Such a solution can always be 
represented as a linear combination of the functions w\^^'^\r), defined by (|2.f 6p - (|2.f 7p . These functions 
have the following asymptotic behaviour as r 0: 



Therefore, the solution vjl'\r) can be written as 

■w\^\r) — cos T]Wi^^'^\r) + sin rjwl^^'^^r), (2.27) 

Note that for special values of SAE parameter, 8 — — {~'^) we have ry = (resp. rj = —). Explicit 

dependence of ?7 on 8, /i, 6 in the general case can be easily found from (|2.2ip and (|2.25p - (|2.26p . Now, 
analogously to the above, consider the following ansatz for the Green function: 



' ' l^^(r) (g) {wl''\r')y for < r' < r < 1. 



T 

for < r < r' < 1, 

(7) {I 
E,l 



This ansatz automatically solves the equation (|2.23p for r ^ r' and satisfies the appropriate boundary 
conditions at the points r = and r = 1. The jump condition at r = r' will be satisfied provided we have 



The last condition trivially holds if one chooses 



f - 
Rr 



A'-^\ = ^ ^. (2.29) 

2 ^/2sin(r7+f) ^ ^ 

Hence the formulas (|2.27p - (j2.29p give the radial Green function for / G (— f , 0). 



2.4 Spectrum 

The spectrum of the full Dirac hamiltonian H^'^^ consists of several parts: 

• a continuous spectrum: > + Ah^ / R^; 

• a finite number of infinitely degenerate Landau levels, given by 



where n — 1,2, .. .Umax < \b\- The allowed eigenvalues of angular momentum correspond to 
-f , -2, ... for 6 > and ?o = 1, 2, . . . for 6 < 0. 
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• a finite number of bound states with finite degeneracy, whose form depends on the sign of the magnetic 
field. Namely, for 6 > one has 



where n = 1, 2, . 



< b — (1 + I/), and for & < we obtain 



^m' + j^[b'~{\b\~n + ^f], 



with n = 1, 2, . . . , n 

& > 0) and lo = -1, 



ax ^ \b\ + ^- The allowed angular momenta are given by Iq — 1, 2, 
-2,...,-<„, (for6<0). 



(for 



• a finite number of non-degenerate bound states, corresponding to the mode with Iq — 0. These energy 
levels are determined as real roots of the equation 



(771 + E)R 



mR 



l+2u 



T{p-b+l)V{^i + b-v) 



tan 



e 



r(-^) 

Note that for ^(8,7^) < we can have a solution of this equation satisfying \E\ < ni. This is in 
constrast with the previous cases, where all energy levels lie in the interval < < + 4b^/R^. 



2.5 Full one- vortex Green function 

Once the Green function of a particular SAE is found, one can also obtain it for any other extension using 
Krein's formula. This fact very much simplifies the analysis of the (S-interaction hamiltonians (see, e. g. [2]), 
since in that case the family of SAEs usually includes free Laplacian/Dirac operator, whose Green function 
can be computed relatively easily (for example, in the planar case one has just to apply Fourier transform). 
The situation with AB hamiltonians is different. Here, the calculation of the resolvent constitutes a non- 
trivial problem even for distinguished values of the extension parameters. 

In the present subsection, we obtain integral representations for the Green function 



G{z,z') 



27r ^ 



e*'^ 1 , 





p — i(io + l)y' 



(2.30) 



of the full hamiltonian H^'^'> for two values of SAE parameter, namely, for 8 = ± — . The outline of the 

calculation is similar to [21] and the reader is referred to this paper for more details. 

We begin by introducing two classes of solutions of the Dirac equation on the disk without AB field: 



/ 



*±(z,( 



±C±e 



9/2 



(l + ze-e)i±A'-''(l-h^e«)±A'+'' 
(l + ze-»)±A'-''(l + ze«)i±A'+'' / 



These functions are delimited by two families of branch cuts in the 0-plane: oo + i{ip + n + 27rZ), Inr + 



and 



■ lnr-|-i((/9-|-7r-t-27rZ), oo + i{ip + Tr + 2TTZ)j , with the arguments of l-l-ze and 1 + ze 

equal to zero on the line ImO = (p. It is also convenient to introduce the "conjugates" of these solutions, 
defined by 

i'±iz,9) ^^±{z < — >z,0< — > -6*). (2.31) 
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The relation L{z)'^±{z,6) — —de'i'±{z,9) allows to construct multivalued radial solutions of the Dirac 
equation with specified monodromy as superpositions of '^±{z, 9). One can check that 

w|^)(z) / e('+V2)«^_(,,^)rf^ = ( ) wl'\r), (2.32) 

JCo{z) \ / 

^(//,±)(^).e/ ^ J^^^^^,iiW^^^,^0)de^2me-^^'l" ^.(m), (2.33) 

where the contour C+(z) (C_(z)) goes counterclockwise around the branch cut ^— oo+i((/3+7r), \nr+i{ip+7r) 

(resp. — Inr + i{ip + tt) , oo + i{tp + 7r)j) , and the contour Co{z) is the line segment joining the branch points 
± Inr + i((y3 + vr). Conjugate solutions are obtained analogously: 

M'\z) f e-('+V2)«^_(,,,)rf, ^ ( ^"'^ 0^ \ (2.34) 

J Co(z) \ / 

w("^±)(z) T / e-('+V2).^^(,,,)rf, = 2,,,-... e-'- ^0^^^^ \ ^^.35) 
'^0^(2:) \ / 

vr 

Let us assume the regularity of the upper component of the Dirac wave function, i. e. ^ ^ ^2' ■ Then 

the Green function (|2.30p can be conveniently expressed in terms of the radial solutions (|2.32p - (|2.35p as 
follows: 

G(z, z') = ^""g"'/' + gi-) /)] , (2.36) 



where 



Q^^\z,z') = for|z|>|z'|, (2.37) 

g(^\z,z') = ^\'''^\^')'»{M'^iz')Y for|z|<|z'|. (2.38) 

Remark. For Q = obtains similar representations, but in this case the summation in fj(^^(z,z') is 

TT 

over ^ ^ — 1. Further calculation is also completely analogous, so we will continue with ® ^ ^'^i ^-"^^ present 

TT 

only the final result for = — at the end of this section. 

2 

Following |21], we substitute into (|237 l) - ((2:38)) instead of w|'^^(2:), w'l'\z'), wl"'^\z), w[^^'^\z') their 
contour integral representations (|2.32p - (|2.35p . After interchanging the order of summation and integration 
the sums over I are reduced to geometric series. For example, in the case \z\ > \z'\ this gives 

gl''\z,z') + gi-\z,z')= j do, J de-2 vi/„(z, e,) ® i>l{z', o^) ""^'Jj^tlf^' /'^ , (2.39) 

Co(z) C+{z')uC-{z') 

where the contours C±{z') satisfy additional constraints: Re(6'i — 6*2) < for all G 6*0(2;), 62 E C-{z') 
and Re(0i — 92) > for all 9i E Ca{z), 92 E C+(z'). After a suitable deformation of integration contours one 
can obtain a representation, which is valid not only for \z\ > \z'\, but for all z, z' such that — ip' ^ ±7t. It 
has the following form (cf. with (2.24)-(2.26) in 27j): 

e-'''('^-'^'+2'^)G(°)(z,z') + A(z,z') for if -if' E (-27r,-7r), 
G(z,z')=<( e-'''('^-'^') G(°)(z,z') + A(z,z') for ip - ip' E {-n , tt) , (2.40) 
e-'^'^iv-'p'-^^) G(°) (z, z') + A(z, z') for - £ (tt, 27r) , 
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with 



in 



Co(z) 



1 - e- 



dOl Ide2 *-(z,0i)®*^(z',02) 



(2.41) 



(2.42) 



Ca(z) lra0:2=v' 



Remark. Starting from (|2.38p , similar results can be obtained. More precisely, one finds again the formula 
([2^ . but with 



G(°)(z,/) = 



An 



dd ■^+{z,6)®^t{z' ,6) 



A(z,z') = -i;2 e-*'^('^-'^') 



1 - e 



8i7r2 



d^l /d02 *+(z,02)®*i(^',ei) 



1 - e^^-^i 



(2.43) 



(2.44) 



Co{z') Im6l2=v 



The proof of equivalence of the representations (|T4T|l . (|2^ and (|T43ll . ((2^ for G(°)(z,z') and A(z,z') 
is left to the reader as an exercise. 

Since for v = Q the second term in (|2.40p vanishes, G'^''^(z,z') coincides with the Green function of the 
Dirac hamiltonian without AB field, which we will denote by H''^\ Using the technique described in the 
Appendix A of [21] . one may compute the integral for G'^'^^{z^ z') in terms of hypergeometric functions. The 
result reads 



G(°)(z,z') 



1 - zz' 
1 - zz' 



-b 



( fl- 



1/2 



1 - zz' 
\l-zz'\ 



Cii(u(z,z') 



ll-zz'l 



Cl2{u{z,z'fj 



C2l{u{z,z')) -[^^) ' C22{U{Z,Z')) 



-1/2 



(2.45) 



where it(z, z') — 



1 - zz' 



and 



an) 



— (1 - m) 2 X 



2nR 



r(i + 2^) 



X 



G+ 2i^i {^l-b+l,^l + h,l + 2^1,1 ~u) 2Fi{fi-b,ti + b,l + 2^1,1- u) 



\ 2Fi{fi-b,fi + b,l + 2fi,l~- u) Cl 2F1 {^1 - b, fi + b+ 1,1 + 2^,1 -u) 

One may also obtain a more explicit expression for A(z, z'): 

-b 



A{z,z') 



smniy 



n 



de 



X 



1 + rr'e'' ^ ^'^ 



1 + rr'e 



Cii(w(r, r',6') 



V 



(1 + r^r'^ + 2rr' cosh 6*) 



1/2 



C2i(i'(r,r',0)) 



l + r2r'2 + 2rr' cosh 6*) 
re~^ + r' 

1 + rr' e^^ 



1/2 



(2.46) 



■e ' 

-1/2 



^'Ci2(^'(r,r',0) 
C22(w(r, r', 



where we have introduced the notation 



v{r, r' ,9) = 



^2 _|_ j,/2 _|_27-r' cosh 6* 
1 + r^r'^ + 2rr' cosh 



(2.47) 
\ 
/ 

(2.48) 
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Remark. In the case © ~ on.*^ finds essentially the same answer (|2.40p . The only difference is that we 

should replace i/ by — 1 in the double integrals entering the definitions (I2.42p and (I2.44p of A(z, z'). This 
finally gives 



A(z,z') = / dd 



1 + rr'e^ 



1 + rr'e^ 
1 + rr' e~ 



1/2 



+ r^r'^ + 2rr' cosh^ 



V 



C2i(i'(r,r',0) 



{l + r^r'^ + 2rr' cosh6l)^/^ 



,1/2 



1 + rr'e 



(2.49) 



Representations (I2.47P and (|2.49p for the vortex-dependent part of the Green function constitute the main 
technical result of this section which will be used later in the construction of Painleve VI transcendents. 



3 Two-point tau function 

3.1 General setting 

In this section, we consider Dirac hamiltonian (|2.3p with two AB vortices located at the points ai,a2 G D. 
Corresponding vector potential has the form 



A 



iBR^ z dz ^ z dz ivi ( dz 



\-\z\ 



dz 



z ^ a\ z — a\ 



iv2 ( dz dz 

2 \ z — a2 z ~ 0,2 



(3.1) 



and it is assumed that — 1 < v\-2 < 0. Let us make a singular gauge transformation H i-^ jji"-^^) — UHU^ 
with 



(C/V^)(z) 



, i^i/2 / \ 1^2/2 

z — ai \ I z ■— a2 ^ 



z ~ ai 



z - 02 



V'(z). 



(3.2) 



It is easy to check that the local action of iJ^"-"^) coincides with that of the free hamiltonian H'^'^^ (i. e., in 
the absence of AB fluxes). However, the functions from the domain of /ft"-") are multivalued: they pick up 
a phase given by e^'^*'^^ when continued around aj {j = 1,2). One should then introduce two branch cuts 
£1, £2 on D as shown in Fig. 1. We do not fix the branches of fractional powers in p.2p . since this will be 
implicitly done later. 













/ r 
























1 



Fig. 1 



Let us isolate the branch cuts in the union S of two open strips S\ and 52 (see Fig. 1) and denote by 
7^(5*) the set of all hamiltonians ijf"^'') satisfying a\ G 5i, 02 G ^2. In other words, the elements of 7^(5*) are 
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parametrized by the positions of AB vortices, their fluxes being fixed. Consider the locahzations of different 
elements of 7^(5') to D\S. Since aU of them are given by the same differential operator, all dependence on 
a is encoded into the spaces of boundary values on dS of the functions from their domains. Several such 
boundary spaces will be considered: 

• a suitably chosen space W — H^/^{dS) of C^-valued functions on dS. 

• a subspace W"^*{a) C W, which is composed of boundary values of functions ifj e H^{S\{£i U ^2)) 
solving the equation (ff^"'"') - E)il; = on S*. Note that W^'"*(a) = W^"'\ai) © W^"'\a2), where 
W™*(aj) {j — 1, 2) is composed of the boundary values on dSj of local solutions to the Dirac equation 
on the strip Sj. 

• similarly, W'^^'^ CW is defined to be the space of boundary values of ^-solutions of (ijf"^") —E)^p = 
on D\S. This subspace clearly does not depend on a. 

• it is also convenient to fix two points €E 6*1, 02 G ^2 and to introduce a reference subspace W^*"*(a"). 

It will be assumed that near each branch point one of the components of the Dirac spinor is regular, and 
the other is square integrable (i. e. we consider four possible types of boundary conditions for H^°'''^^). 
With this choice of the domain, the spaces VF*"*(a) and VF*^^* can be shown to be transverse in W for 

E e C\|^(— 00,— m] U [m,oo)^. Then, given any subspace V C one can define the projection P{a) : 
V ^ W"'\a) of V on W"'\a) along W''^. 

Proposition 3.1 Let f ^ W and consider the function fint G H^{S\{£i U £2)) defined by 

f^nt{z) - - / &'^'''\z, z')-^^ {a_ f{z') dz' + a+ f{z') d-z'} , (3.3) 
JdS ^~ \z \ 





{ 






f 




-+ = l 


V 


; • 


- = l 




I) 



where z G S\{£i\j£2), dS is oriented counterclockwise and &'^''^\z, z') denotes the Green function of H^°'''^\ 
Then the boundary value of f^t on dS coincides with the projection P{a)f £ W"^*{a). 

M First remark that fint satisfies Dirac equation on S. Therefore, we may write 

f^nt{z) = / G('^''')(Z,^') (^-■•'^ - e) f,nt{z') df,,, . 

Jd\s ^ ' 

Integrating once by parts and using Stokes theorem, one obtains 

f^nt (Z)^ - [ {Z, Z') f^nt {z') dz' + a+ f,nt {z') dz'} , (3.4) 

JdS ^~ \Z \ 

and thus the map (|3.3p is indeed a projection on T^*"*(a). Using similar arguments, one can show that its 
kernel coincides with . □ 

Definition 3.2 Let fj e H^/^{dSj) (j = 1,2) and consider the function fint.j & H^i^jX^j) defined by 

hntA^) = - I G^'^^'''^\z,z')-^^{a^f,{z')dz' + a+f,{z')d-z'}, (3.5) 

where z G Sj\£j and G'^°'^'^^\z, z') is the Green function of the Dirac hamiltonian on the disk with only one 

branch point a j . Passing to boundary values of fint, j ondSj, one obtains a projection Pj{aj) : H^^^{dSj) *■ 

W™*(aj). We will denote by F{a) the direct sum of such one-point projections: 

F{a) =^ Fi(ai) ^2(^2) : W M^*"*(a). 
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Definition 3.3 r-function of H^"^''^^ is defined as follows: 

T(a,a°) = det f [Pi(ai) Paloa)! W'"\a°)] = detw^^Ha") {P-\a)F{a)) . (3.6) 

y L J H/'"* (a") / 

Remark. For this definition to make sense, the restriction of the map P^^{a)F{a) to M^™*(a°) should be 
a trace class perturbation of the identity on W^™*(a°). This can be shown analogously to Proposition 4.2 of 
[27] . Factorization formulas for the derivatives of one-point Green functions {z, z'), needed for the 

proof, are presented in Subsection 13.4.31 



We are interested in explicit calculation of the r-function p.6p . In order to do this, it is convenient to 
introduce yet other boundary maps. 

Definition 3.4 Let j be a smooth curve dividing Poincare disk into two disconnected parts Di and D2 as 
shown in Fig. 2. For any f G H^^'^{'-f) define the functions 

f±{z) ^±1^ G^°\z, W- fiA dz' + C7+ f{z') d-z'} , (3.7) 

where z £ D2 for f+{z), z £ D\ for f-{z) and G^^^{z,z') is the Green function of the Dirac hamiltonian 
without AB field, explicitly given by \2.45^ - {2.46^ . With some abuse of notation, the boundary values of 
f±{z) on 7 will also be denoted by f±{z). 







\ 






1 






1 



Fig. 2 



Note that f± have the properties of projections: 

(/±)±=/±, (/±)^ = 0, 

and, in addition, /+ -I- /_ = /. One can thus uniquely write any / e H^/'^{'^) as a sum of two functions, 
where the first term (/+) may be continued from 7 to D2 as a solution of the Dirac equation without AB 
field, and the second one (/_) is the boundary value of a solution on Di. 

Let us represent any function -i})'^^^ e H^^^{dSj) {j — 1,2) in the following form: 

',0) - 









r ) 







where the indices L and R correspond to the left and right boundary of the strip Sj . 

Proposition 3.5 If^^J^ is the boundary value of a solution of the Dirac equation on Sj without AB vortex, 
then its components satisfy the relation 



^U)^) \u{dSf\dSf) 



(3.9) 
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where the integral operators d) {dSj'\dS^^ , uj {dS^\dS^^ are defined as follows: 

{dSf\dSf) V'i^l) (z) = - [ G(°)(z, z')-^^ |a_ ^«_(z') dz' + a+ i^ll{z') dA , (3.10) 

can be continued to a solution on Sj, then, similarly to (|3.4p . we have 

i,^^\z) = - [ G(°)(z,z')-^^|a_^(j)(z')dz' + t7+^(^)(z')dz'| = 
Jas, 1 - FT J 

iR 



where the second equality follows from the fact that it is not possible to construct a solution on the whole 
disk D, belonging to the domain of the free Dirac hamiltonian H'-^\ Passing in p.l3p to boundary values 
and taking the projections, we obtain the result p.9p . □ 

Proposition 3.6 If tjA^^ is the boundary value of a solution of the Dirac equation on Sj\£j with only one 
branching point Oj, then its components satisfy the relation 



(3.14) 



(/35,K)^tL)(z) = -y^^^G("-^^)(z,zOY^^{'^-^i^?-(^')rf^' + '^+V4-(^')dz'}, (3.16) 
(7s,K)^i^;V) W = -| ^G('^-''^-H^,zOY^^{<T-^i^,^+(zOrf^' + '^+V^l(^')rf^'}, (3.17) 
('55,K)V4-)W = -/ ^A('^-^^-)(z,zOY^^{'^-^l^?-(^')d^' + ^+^i^l(^')rf^'}- (3.18) 



where the integral operators asj{aj), Psjio,j) ISjio-j)? ^Sjio-j) defined as follows: 

(«s,(a,)V'^;V)(^) = -^^^A('^-''^)(z,zOY:^{cT-^i';V(^orf^'+^+V'i';i(^orf^^ (3.15) 

iR 

IdSf 

iR 

ij"-^ {z, z: ) — 

IdS'r 

iR 

'asf 

Here, we denote M''^'''^\z, z') = G'^'^^'^^^z, z') - G^^^z, z'). 
■ The proof is completely analogous to the previous one. □ 

We now choose + cind ipR^- to be the "coordinates" in W™*(aj) (j = 1,2). It is not difficult to 
check that in these coordinates the map F{a) : l/l^™*(a°) iy™*(a) is given by the identity operator. In 
order to find the representation of P{a)~^, one should be able to decompose any g S Vl^™*(a) as g ^ f — h, 
with / £ 14^™*(a'') and h € W^^^ . More precisely, we are interested in the relation between / and g. This 
calculation is very similar to the proof of the Theorem 3.1 in [27]: 

• First note that ft,^^' can be continued as a solution to the left of dS\^ and, analogously, can be 
continued to the right of 95*^. Then one has — /i^L — and, therefore, 

^^^V = /i!l ^-=f^^- (3.19) 
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Next recall that the boundary values and h^^' are not independent. They are related by the 
formulas 

;,(2) 



,(2) 



wi2 = u; (dS^ldS^) , UJ21 
which follow from the Proposition 13. 51 Then one finds 




(2) 



/•(I) 

f(2) 



r, f(2) 
r, f(i) 



• Finally, according to the Proposition 13.61 one has 



Af) 

J L- 

A3) 
J E 



9^. 



7s,(«°) <5s,(a°) 




(3.20) 
(3.21) 

(3.22) 
(3.23) 



If we now find g^]^\, f^^, ffi^L, /l^L from (|3.22p - p.23p and substitute the corresponding expressions into 
p.20p - p.2ip . two more relations between the coordinates of / and g can be obtained. Together with p.l9p . 
they may be written as follows: 



/ 



V 



1 



-'^2i7Si(ai) 


1 



-'^2i7Si(ai) 






1 

-W2i5si(ai) 




1 

-W2ii5si(ai) 






-cJi2ds2(a2) 
1 




-Lbi2as^{al) 
1 






-^12/352(02) 


1 



-'^i2/3s2(a^) 




9 b'- 
9Z 
V 9't ) 



(3.24) 



f(2) 
J L.+ 



J 



or, in a more compact form, 

(1 + Af(a))g= (1 + M(a"))/. (3.25) 

Therefore, the r-function is given by 

T(a,a°) = detvi.-*(aO) (l + M(a°))"^ (1 + Af(a)) . (3.26) 

Remark. One-point Green function &'^^'^^\z^z') is related to the one-vortex Green function G{z,z'), 
calculated in the previous section, by a simple unitary (-1- singular gauge) transformation. For example, one 
has 

^ia,=o.^,=^)^^^^>s^ = e*''('^-'^')A(z,z')- (3.27) 
Thus the action of all operators in p.24p is known. However, in order to obtain from p.26p an explicit 
formula for the r-function, one still needs to introduce some coordinates in the infinite-dimensional spaces 
(^H^/^{dSlf)J ^. It turns out, though, (see Theorem 6.3 in [29] and Subsection 13.4.31 of the present paper) 

that the logarithmic derivative of the r-function (|3.6p does not depend on localization (i. e. on the choice of 
5) and on a°. Moreover, since r(a°,a'') = 1, one should have 

r{a) 



r(a, a 



Therefore, one may fix the coordinates in (H^^^{dSi'^)] and calculate the r-function for a convenient 



choice of S. These problems are addressed in the following subsections. 



16 



3.2 Dirac equation on the Poincare strip 

Coordinate change z = tanh^ maps the Poincare disk onto the strip U — |c ■ l^yl < ^| i'^ the complex 
^-plane. In these coordinates, the Poincare metric is given by 



= g^^d^di = 



didi 



cosh^ - e) ' 

and the vector potential of the uniform magnetic field B can be chosen in the form 



^^^^4--26tan24dC.- 



Corresponding Dirac hamiltonian explicitly given by the formula (|2.3p with 

K = |2 cos 2^5^95 - i(l + 26) sin2^y|, 
K* = --^ |2cos2fyaj- + i(l - 26)sin2^y|. 
It is related to the free Dirac hamiltonian on the disk by a unitary transformation 



(0 = UosiO H^;^i^ ^ tanhO Ul^iO 



(0) 



^ / cosh A ' 2' 
\ cosh f / 



/ cosh^ 
\ cosh ^ 



(3.28) 

(3.29) 
(3.30) 

(3.31) 
(3.32) 



The main advantage of the gauge (|3.28p is that H^^^^^ commutes with the ^a:-niomentum operator = —id^^ ■ 
The eigenspace of Px, characterized by the momentum p g K, is composed of the spinors of the form 
g(p, ^y)e'P^=" . Being restricted to this eigenspace, the hamiltonian ^g^^^'ip acts as follows: 



g(jP,iy)^ Hpg{p,^y), Hp = R M 



where the operators Kp and K* are given by 



Kl = -i 



cos 2^y {d^^ - p) + (1 + 26) sin 2^y 
cos2Cy {d^^ +p) + (1 - 26) sin2^y 



Let us consider the partial Dirac equation 

{Hp-E)gip,^y) = 0, 



(3.33) 

(3.34) 
(3.35) 

(3.36) 



where it is assumed that E is real and \E\ < m. Two linearly independent solutions of p.36p can be chosen 
in the following way: 



$(±)(p,^,) - (2cos2e,)— Vxb) X 



(3.37) 



\ 



I 
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with 

r{^,-b + i)r(/i + 6 + i)r + 1 - f ) r + 1 + f ) 

47r(r(l + 2Ai)f 

The "~" in 2F1 indicates that the hypergeometric function is defined not on its principal branch, but on the 
cut plane C\(— oo,l] with Hm 2Fi{a,b,c, z) = 1, as in the formula (3.5) of [3. Notice that the first 

solution, (p, S^y ), satisfies the condition of square integrability on the upper edge of the strip U {£_y ~ —) 
and the second one, <1>(^^ (p, ^j,), does so on the lower edge (i^j^ = 

{p,^y ^ ±J) = (2cos2C,)'^ VX(P) ( ) +0 ((cos2^,)'^) . (3.38) 



These two solutions verify simple symmetry relations: 



<^>(+){p,Q = a,^^~){-p,-(y), $(±)(p,C,) = a,$(±)(p,^^). (3.39) 

It is also worthwhile to give a formula for the determinant of the fundamental matrix built from ^^^^ (p, £^y) 
and $'^^)(p, ^j,). Using transformation formulas for hypergeometric functions, one obtains 

det($W(p,ej,),$(")(p,Cs,)) =-icos2iy. (3.40) 

Green ftmction Gs.pi^y, £,y) of the partial hamiltonian Hp satisfies the equation 

[HpiO - E) GEAivA) = ^{iy - Q I2. (3.41) 

Analogously to Subsection 12.31 consider the ansatz 



^ /c^.p$(-Hp,e.)®(<i'(+)(P,e;))' for -f <e,<e;<f, 

\Ce,P $( + )(p,Cy) ® [^^^-\vXy)) for -\<^y<iy<\- 

It solves (I3.41[) for 7^ ^y and satisfies the appropriate boundary conditions as ~^ • Required 

singular behaviour at i^y = is equivalent to the condition 

cos 1^ 

GE^iv + 0, iy) - GE.Aiy - 0, ^y) - z — ^ a,. (3.43) 

Substituting (|3.42p into the last relation and using symmetry properties (|3.39p . one may show that (|3.43p 
holds true if we choose Ce,p ~ R^^- 

3.3 Fredholm determinant representations 
3.3.1 Boundary projections revisited 

Let us fix a hue C^(a) = \^ 'EUlS.y = r ^^id consider a C^-valued function gAo}(^x) G i7^/^(£,(o) ), 
represented by its Fourier decomposition 

/oo 
dp g{p,il°^)e'P^'^. (3.44) 
-00 

Next we introduce two operators, Q±{£,y^'^), whose action on Fourier transform is given by a matrix multi- 
plication 

o±(?r)5(p,4"^) = Q±(P,4°^)5(P,4"^), 
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with 



(3.45) 



Symmetry properties p.39|) and the formula (j3.40|) imply the following result: 
Proposition 3.7 The operators Q±ij>,S}y^) satisfy the relations 



(3.46) 



Thus one has a decomposition H^^^C^w) = -ff+(4°^) © ^-(4°^)' ^l^^re i?±(C^°^) - Q±iC'y'^)H^^^iC^m). 
We remark at once that -ff±(^^"-') will play the role of the subspaces (_ff^/^(9S'f 2^) ) from Subsection 13. 1[ 



since they are composed of the boundary values of functions satisfying the Dirac equation (jT^^^-^p — tp = 

and the condition of square integrability in the region > (in the case of ff+) or < (for H^). 
Using the formulas ()3.45p - (|3.46|) . one may write Fourier transform of gAo}{^x) as follows: 

g{p, ) = 9+{p, d°))$(+) (p, d")) + 5- {p, C')<f>(-' {p, e^), 



where 



9±iP,i 



(o)^ _ 



cos2C^''^ 



(3.47) 



It is convenient to think of ,g± (p, ^^"^ ) as coordinates of 5,(0) in i?±(Cy°''). It should be noted that g±ip,^y^^) 



are ordinary functions, in contrast to p.Sp . where a similar notation was used for 2-columns. In such 
coordinates, the Propositions 13.51 and 13.61 have the following form: 





71/4 






-^-y 


z 


—e 








e(L) 
^ y 




-71/4 



Fig. 3 



Proposition 3.8 Let us consider a strip 5 = |^ G | < S^y < ^y^^^- Suppose that a -valued function 
can be continued to S as a solution of the Dirac equation i^H^^l^^ — ^ = 0. Then one has 



i^R.Ap) ) V 1 y V 



(3.48) 



Proposition 3.9 Let us now assume that the strip S contains one branching point ^° and introduce a 
horizontal branch cut £ £ S as shown in Fig. 3. Suppose that tp S H^^^{dS) is the boundary value of a 
multivalued solution of the Dirac equation on S\i, characterized by the monodromy e^'^"' at the point 
Then 

^lAp) \^( "5(r) /35(e) \ ( AAp) 
i^uAp) J V 75(e) 4(e) J V i^nAp) 
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where 





){p) 


= R 


— oo 




-{q) dq, 


(3.49) 




)ip) 


= R 


— OO 


''Hp, q) ^r,- 


-{q) dq, 


(3.50) 


(75(r)v^L,+ 


)ip) 


= R 


— oo 


"Hp-, q) i-'L,^ 


(g) dq, 


(3.51) 




)ip) 


= R 


— OO 


"Hp, q) ■4'R,- 


.{q) dq. 


(3.52) 



and 



Af^''\p,q)^:, ^7 (3-53) 

271 COS 2^y COS 

1 

27r cos 2^y cos 2^(, 

d^dC e-^«^+'< ($(±)(p,e,))V. Gii;,^^(e,e') ,,,,,, 'y.^^^\q,Cy 



— 00 —00 

Gf '"^^p,?) = — X (3.54) 



— CO —00 



Here, G^strip\^^^') d^'^otes the Green function of the Dirac hamiltonian on the strip with one branching 
point e, and A^;;^(^,C) = G''£i;\tC) - G%{C,C)- 

■ We illustrate the idea by deriving (|3.49p . On the strip, the formula (|3.15p is replaced by 

(a5(r)V'L,+)(e) = — ^ r Aii:^;\^,e)<7^i'LAadc, 

^ ' cos2Q "'-00 

with ^y — ^y— ^y"^ . Then after Fourier transform one obtains 

(a5(r)^L.+)(p)-<i'^"'(p,4^^)= (3.55) 

(00 00 '\ 
^ / A£p'')(e,e')e-^^«^+^^«4a.<i>(+n9,4'^^)^L.+ (9). 



cos2ef ^ 



— 00 —00 



The columns of the expression in curly brackets satisfy partial Dirac equation (iJp(^) — £')V' = and, in 
addition, for ^y < they are square integrable as — — . Hence they are both proportional to {p, S,y). 

Similarly, the rows of this expression are proportional to ($(^) ((7, ^y))^. Therefore, for S,y,£,y < one has 



00 00 

1 



2n 



d^.dC Alf°p'')(^,e')e-'^«==+^««i = A^f^''\p,q) ^^-^P.^y) ® (<f ^"^9, C;)) \ (3-56) 



-00 —00 



and it is a simple matter to verify that the coefhcient A*^ '''\p,q) is indeed given by p.53p . Substituting 
([336)1 into (lO^)) . we obtain ((09| . □ 
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Finally, suppose that the Poincare strip contains two branching points, and ^"2, such that tanh^^J = 
(j = 1, 2) and <^°i < ^"^^ Using p.3ip and going through the definitions of Subsection l3.11 one finds that 



the two-point tau function r(a) can be written as a Fredholm determinant: 

r(a)=det(l-d(e)^(e)), 



(3.57) 



The kernels p.49p and (|3.52p of the integral operators and S{£^°-'-) do not depend on the choice of the 

strips Si and hence the corresponding indices will be omitted from now on. 

Remark. It turns out that the tau function depends only on the geodesic distance d{ai,a2) between the 
points ai and a2 (see Appendix). Therefore, to simphfy (|3.57p . we may choose = 0, ~ Ig + iO, 



I. 



(01,02) 

R 



In addition, the invariance of H^^^-^^ with respect to ^a:-translations implies that 



Thus it remains to compute only the quantities A'£''^\p,q), henceforth referred to as "form factors". This 
problem is solved in the next subsection. 



3.3.2 Asymptotics of A(^,$') 

Let us first compute A^^''^\p,q). The idea is to use instead of (|3.53p the equivalent formula (|3.56p . Since 
it is valid for all ^y,Cy < 0, one may obtain A^^''^^ (j),q) by exploring the asymptotics of both sides of this 
relation as ^y,^'y — — ■ The asymptotics of the RHS of p.56p is determined by the formula (|3.38p . which 
gives 

-I c\ 

(3.58) 

TT 

To find the asymptotics of the LHS, we will use the results of Section [5J Set for definiteness ® = ~ 2 ' 
The formula (|2.47p then implies that to leading order 



leading behaviour a(0-'')/ \ I 1 \ t \ Ia nt nt'\- 

of the RHS of & as iy,Cy - -f ^ '^^ Vx{p)x{q) (4 cos 2^, cos 2^,) 



A(z,z') 



sinvTiy r(/i - + l)r(/i + 5 + 1) / C 



r,r' — >-l,\ip—ip'\<C7r 



r(i + 2^) 



-tip 



(3.59) 



X [(l-r^)(l-r 



/2\1 2 



de 



e 2 



e'+^(^-^'^ + l (2 + 2cosh0)^ 



In order to obtain from (|3.59p the asymptotics of Ag°^|^p (f , <^') as ^j^, — — , one should take into account 
the following comments: 











■ll 


i 1 — 

11 









nIA 




I 










-ji/4 



Fig. 4 



21 



• Dirac hamiltonian on the disk with one AB vortex is related to the free hamiltonian on the disk with 
one branching point by a singular gauge transformation, see p.27p . The choice of the branch cut 
represented on the Fig. 4 is equivalent to saying that (p and ^p' vary between — tt and tt. 

• Free hamiltonians on the disk and on the strip are related by the unitary transformation (j3.32p . Also 
notice that 

' e»^(^+f) \ 



Uds{^ arctanhz) 







ip e (-TT, 0). 



The integrals over and in (j3.59p can be written as integrals over ip and p' using the formulas 

1 - dip 



ctg(-|), 2 cos 24 



sin(— (y9) 



f 2sin(-93) ' 



where, as above, p} G (— tt, 0). 



Once the leading behaviour of A^l'jlj'^p (^, is found, one should substitute it into (|3.56p . Then, comparing 
the LHS of (|3.56p with p.58p . we obtain a triple integral representation for the form factors: 

R A^^'"^ {p, q) = Vpip)p{q) {p, 9) , 

where 

22^r(i + 2^) 



p{p) 
^u(p, q) 



r(;. + i + f)r(M + i-f)' 

00 7r/2 7r/2 

dO I dx I dy 



smTTZ^ 



-00 



(^^e-2t{x-y) ^ 1) (2 + 2cosh6')~ 

X (sinx)^+"f ~^(cosx)''^'^^^(sin?;)^^^^5 (cosy)^ 
One may deduce from p.60p - p.62p the following symmetry properties: 



(3.60) 
(3.61) 
(3.62) 



{p,q)^A^^'-\p,q)=A'^^''>{q,p) 



(3.63) 



Another useful representation for A^^''^^(j}, q) can be found by decomposing the integral over 9 in (|3.62p 

into two pieces, and , and replacing (^e^-'^^i'^-v) + ^ by the appropriate geometric series in each 
piece. Then the integrals over x, y and 9 can be computed in terms of hypergeometric functions and we 
obtain 

1/2 



RA!~^-''\p,q)^2'^^-^ 



r(M + i + f)r(A. + i-f)r(^ + i + H)r(^ + i-f) 



' (^_l)n-lg-7r(p+9)/4 

/ ; ; 2F1 

^ u+l + iy-b + n 

^ n—l ^ 



2F1 
2F1 



^ (-l)"e^'^''+'?)/'* 
r^— — 2F1 



a - v + h + n 

71— ' 



2F1 
2F1 



r(i + 2^) 



Ai + l + i^-6 + n,^+i + y,l + 2^,2-i0 
1 iq 

/i + l + i/-6 + n,/i+--y,l + 2^,2 + i0 



^ + l + z/-6 + n, l + 2/i,/i + 2 + 1^-6 + + 



^-t/ + 6 + n,^+i + y,l + 2/i, 2 + iO 
1 iq 



II — v + b + n,l + 2/i, fi — ly + b + n + l, — 1 



X (3.64) 
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It is worthwhile to note that both series in (|3.64|1 rapidly converge, since for large n their n-th terms tend 
to as at least n^^^^*". 

Remark. We have obtained the representations (|3.60p - (|3.62[) and (|3.64[) assuming that p,q G M.. However, 
they also give an analytic continuation of A^''^-'(p, q) to the strips |Imp|, \lmq\ < 1 + 2/i. The extension to 
the whole complex p- and g-plane may also be constructed: in particular, it can be shown that !F^{p, q) has 
only simple poles at p = ±z(l + 2/i + 2ni), q ~ ±i(l + 2/i + 2?i2), where n\^i = 0, 1, 2, . . . This singularity 
structure will be used later in the study of the long-distance behaviour of the r-function. 

The calculation of t^^'^\p^q) is quite similar, the result being simply t^^"'^\p^q) — A'"^''^\—p,—q). 

TT 

Finally, for another distinguished value of the SAE parameter, ® = representation (|2.49p for the 

corresponding Green function implies that 

RA^i'''\p,q) = ^pip)p{q)T,-i{TP,Tq). (3.65) 

Thus we can now compute the r-function of the Dirac hamiltonian for all four types of boundary conditions. 

Example. Let us choose = 0, = Is + iO, as in the remark at the end of Subsection 13.3.11 Recall 
that Is — arctanhy^ denotes rescaled geodesic distance between the points ai = and 02 = s. We also 
assume that the upper component of the functions from the domain of the Dirac hamiltonian is regular as 
^ _> ^ ^<»2 ^ Then the r-function p.57p can be written as follows: 

r(s) =det(l-L,,,,L:,^_,), (3.66) 
where the kernels of the integral operators -/ji/,s, L'^^ s 

L^AP^q) = e'(P-9)/./2 ^p(p)p(q)^^(p,5), (3.67) 

l'vAp^i) = L^^s{-p,-q), (3.68) 

and p, g S R. The functions p{p), T^{p, q) are given by (|3.6ip . (|3.62|) . 
3.4 Relation to Painleve VI 

Let us briefly describe the relation of the present paper to the PBT work [29]. Recall that the hamiltonian 
H'^^'^ of a Dirac particle in the absence of the AB fluxes is given by the formula (|2.3|) with 

K = i?-i[2(l-|z|2)a, + (l + 26)z], 
K* = [2(1- |znaj + (l-26)z] . 

Consider the operator 

It is straightforward to check that A coincides with the operator m — studied by PBT (see, e. g. the 
formulas (1.14)-(1.16) in ^) if we identify 

mpBT = \/ — E"^ , A:pBT = — 

In the presence of branch points, one should only replace 7?^'^^ in the definition of A by the operator H^°-''^\ 
introduced in the beginning of this section (recall that is obtained from the Dirac hamiltonian with AB 

field by a singular gauge transformation). Thus there is a unique correspondence between the multivalued 
solutions of the Dirac equation considered in and the solutions of [H'-"-'"') — E)ip = 0. Using this 
correspondence, we now reformulate the key steps of the PBT analysis in the context of the present work. 
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3.4.1 Symmetries and elementary solutions 

The hamiltonian H^'^^ transforms covariantly under the action of the isometry group of the Poincare disk. 
In particular, if iji^z) satisfies the equation [H'^'^^ — E)iIj — 0, then for any g — [ ^ ) G SU{1, 1) the 



function tpg{z) defined by 



^,{z)=\ ^(^''^^ ^ ]^{z^-.{z)), = (3.69) 

v{g,z) 2 1" a-fjz 



is a solution of the same equation. A basis in the Lie algebra of the corresponding complexified infinitesimal 
symmetries can be chosen in the following way: 

M+ = -z'^d, +dg + bz~ ^a^, 
M_ d,-z^dg-bz+^a,, 
M3 = zdz — zdz ^^+2"^^ = L — b. 
These generators satisfy sl{2) commutation relations 

[Af3, M±] = ±M± , [Af+, M_] = 2M3, 

and, therefore, one may introduce two families of solutions {wi}, {w;*} of the Dirac equation on D with one 
branch point at z = 0, on which the symmetries act as follows: 

M3W1 = il-b)wi, M3WI = -{l + b)wl 

M+wi = m+{l,b)wi+i, M+w* = w*_^, (3.70) 

M^wi — wi^i, M^Wi — m^{l,b)w'^j^-y. 

These relations fix {w/}, {w;*} up to overall normalization. Since {w/}, {wj*} diagonalize both the hamiltonian 
and angular momentum, they are related to the radial wave functions introduced in Section 2. It is convenient 
to choose 



( ^ r(^ - 6 + 1) (77,+), X 
"'("^ ^ -^ r(M-5 + ? + i/2) "'-V2 (^)' 



e 



r(^ + & + i) 



"^'^"^ ^ ^ r(^ + b + / + i/2) "-'-v^^"^' 

and it then follows that m±{l,b) — /i^ — (6 =F Z =F 1/2)^- Using (|3.69p . we also introduce the elementary 
solutions with one branch point at z — a: 

wi{z, a) = V{T[a\, z)wi{T[-a\z), w*i{z, a) = V{T[a\,z)w*{T[-a\z), 

where 

/ l-2b _ l + 2b \ 

z — a ,,,^r IS ,. / f 1 — az\ ^ /l — az" 



T[— a]z = V{T[a], z) = diag 



1 — az \\1 — az J \1 ^ az 

3.4.2 Local expansions and deformation equations 

Now consider multivalued solutions of the Dirac equation, which are branched at two points 01,02 G D 
with fixed monodromies e^'^''^^'^. In a sufficiently small finite neighbourhood of each branch point any such 
solution can be represented by an expansion of the form 

neZ+1/2 
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It is convenient to introduce instead of ia,- a new parameter i>j (j ~ 1,2), which can be equal to either i^j or 
+ 1. We assume that —1 < i'i^2 < and thus < pi,2| < 1- Consider the response functions Wj{z,i>) 
and W*{z, v) {j = 1, 2), which satisfy the foUowing conditions: 

• Wjiz^v) and W*{z,i>) are multivalued solutions of the Dirac equation with the above monodromy 
which are square integrable (with the measure dfj,) as \z\ 1. 

• Wj{z, £>) and W*{z, v) have local expansions of the form 

n>0 

W*{z,v)[ak] = 5jk^*_^/2-Oki^-.(^k) + ^[cn^jUJn+u^{z,ak) ^ d^nji^l_^^{z,ak)] (3.72) 

n>0 

where fc = 1, 2 and n e Z + 1/2. 

It turns out that for real values of E such that \E\ < m these requirements fix Wj{z,i>) and W*{z,v) 
uniquely. Therefore, the expansions p.7ip - p.72p can be thought of as defining the coefRcients a^^-, b^^, 
c^ j and d^ j as functions of a and i). 

The lowest order coefficients satisfy a set of deformation equations in a (see Theorem 5.0 in [29]). If we 
introduce the 2x2 matrices with the elements Ajk — cijdjk, Ajk — djSjk, Ajk — S'jSjk and also 









e = 


A - 61 + [ai, A], 


(bi),fc 




KH, 


/- 


AbiA - bi, 


(Cl)jfe 




KH, 




ci — AciA, 


(di),fc 


= 4/2.,/(l- 


kP), 


= 


A-W-[di,A], 



(3.73) 



these equations are given by 



where 



de = fG + Fg+[E,e], (3.74) 

df = eF + Fh + Ef - fH, (3.75) 

= hG + Ge + Hg-gE, (3.76) 

d/i = gF + Gf +[H,h], (3.77) 



E = (A-(6 + l/2)l)^— ^^ + [dA,ai], 

F = dAhiA + AhidA, 

G = d^CiA + ^CidA, 

H = -(A-(fc-l/2)l)^— ^^ + [dA,di 



One also has symmetry relations 



ef~fh = ge-hg = 0, (3.78) 

e^-f9^h^~9f = A^'l, (3.79) 

[esin7rA(l- = hsinnA {1 - \A\^), (3.80) 

[/sin^A(l-|An]^ = /sin7rA(l-|An, (3.81) 

[gsiuTrA (1- lAp)]^ = g simrA {1 - \A\^). (3.82) 
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In addition, the diagonal elements a{/2 j' '^1/2 j be expressed as follows: 

«i/2,j = - 1/2, b) ^3^°'\i2,k + H ^J'="fe4/2,fc' (3-83) 

fe^J fc 

^i/2j = V2,^)-Xl'^Jfe4/2,fc-Hfj'=«fe^i/2,fc- (3-84) 

fc#J fc 

It is known that the system (|3.74p - (l3.77p combined with the relations (I3.78p - (j3.82|) can be integrated in 
terms of a Painleve VI transcendent. In particular, if we choose a\ = 0, a2 = \fs and set 

^ = ^, (3.85) 
/11/22 1 - s 

then for v\^i > the function wis) satisfies Painleve VI equation (|l.ip with parameters (|1.3p (for the details 
of the proof, see [H]). 

3.4.3 Tau function 

The link between the deformation equations and the tau function considered above is provided by a formula 
for the derivative of the Green function of 

(l-|a,f)9„^.G("-)(z,z') = ^ , W,[zrv)®W;{z',v)\ (3.86) 

{l-\a,\^)d-aP^^^^\z,z') = ^„ -L,,. W;{z,v)®W,{z\v)'^. (3.87) 

Remarkable factorized form of these expressions follows from the fact that Green function (z, z') inverts 

the operator — E, and from the analysis of the local expansions of G^"'"-* (z, z') near the branch points. 

Numerical factors in ()3.86|) - p.87p may be determined using a variant of Stokes theorem calculations from 
Subsection 13.11 Different boundary conditions for TJ^"'''^ are encoded into the choice of {vj}: the value 
Vj = Vj + 1 (or i^j — Vj) corresponds to the functions whose upper (resp. lower) component is regular at aj. 

Using (|3.86P " (|3.87p one can show (see Theorem 6.3 in j29j) that the logarithmic derivative of the r- 
function (|3.6p may be written in terms of the lowest order expansion coefficients of Wj{z, v) and W*(z, i>): 

d\nT{a,a°) = X! 1 _ |a.|2 {°i/2,j^°J + g^i/2j'^^j}- (3-88) 

Specializing this formula to the case oi — Q, a2 = ^/s and 1)1.2 > and using (I3.83p - (|3.84p . PBT have 
obtained the relation (II. 2p between the corresponding r-function and the Painleve VI transcendent defined 
by (|3.85p . Analogous result holds for arbitrary branch point positions, as the r-function actually depends 
only on the geodesic distance between ai and 02. The proof of the last statement is given in the Appendix. 



3.5 Long distance asymptotics 

Asymptotics of the PBT r-function as /s ^ 00 (i. e. as s ^ 1) can be obtained by expanding the Fredholm 
determinant ()3.66p : 

lnr(.) = -y _ (3 89) 

n=l 

Leading order behaviour is determined by the first term in p.89p . 



00 00 



Tr(i,„,L:,^,,) - J dp J dq p{p)p{q)T,,{p,q)T,,i-p,-q)e'^P'^^'^. (3.90) 



-OC' —OC 
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The integrand in (|3.90p has particularly simple analytic properties in the complex p- and q-plane. It has 
already been noted above that Tu^^Pil) ^i^d Ty^{—p,—q) have simple poles aX p — ±i(l + 2/i + 2ni), 
q = ±i(l + 2/i + 2n2), where ni_2 = 0, 1, 2, ... At the same points, the functions p(p) and /9(g) have simple 
zeros and thus the whole integrand has simple poles. We can then close the contours of integration over p 
and q in the upper and lower half-plane, respectively (recall that Is > 0). Summing over the residues, one 
finds a long-distance expansion of the integral p.90p . The leading term of this expansion corresponds to the 
poles at p = i{l + 2/i) and q = — i(l -I- 2p) and yields 



Tr {Ly^^sL'^,^ g) 



YesTy^(p, q) 



vesTy^[p,q) 



p = i(l+2tj,) 
q=~i{l+2ti.) 



-2(l+2p)h 



o 



(3.91) 

The residues in the last formula can be extracted from, e. g., the representation (|3.64p . After somewhat 
cumbersome computation one finds 



2sin7ri/ r(/i + 2 + v - b)T{fi -ly + b) 



p = i(l+2M) 
q=-i(l+2fi) 



Combining ((3^ - ((3^ and using that 4e" 



r(2 + 2n) 
1, we finally obtain the asymptotics 



(3.92) 



t{s) ~ 1 - Ar{l - s)l+2^ + O ((1 - 5)2+2^") 

with Ar given by (|1.5p . This finishes the proof of Theorem 11.11 



1, 



4 PVI PV: flat space limit 

In the present section, the analogs of the above results in the limit of flat space are established. Since 
the self-adjoint extensions and the spectrum of the corresponding Dirac hamiltonian have already been 
discussed in [12] (see also [6] for the zero- field case), we will not dwell much on this point. One- vortex Green 
function was also computed in a closed form by Gavrilov et al [14'. However, the representation found in 
p"4j is inconvenient for our purposes, so below we obtain another formula, in which the vortex-dcpcndcnt 
contribution to the resolvent is manifestly separated from the "free" part. This formula enables us to find 
Fredholm determinant representations for the two-point tau function of the Dirac hamiltonian on the plane, 
which turns out to be related to a class of Painleve V transcendents. 

4.1 One- vortex Green function 

Quantum motion of a Dirac particle on the plane in the presence of an external magnetic field is described 
by the hamiltonian 

As above, we will consider a vector potential describing the superposition of a uniform magnetic field B and 
of the field of an AB fiux <& = 27rj/ (—1 < ly < 0), situated at the origin. Analogously to (|2.6[) - (|2.7[) . one can 
choose 

A^^'> = --B{zdz- zdz) ^ —dip, (4.1) 
.(,A ( dz dz\ , , , 

and then the hamiltonian becomes 

H = 



Bz V 




(4.3) 
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Radial hamiltonians -ffzo+y (^o G corresponding to different angular momentum eigenvalues (equal to 
Zo + -), are defined by 

^ m dr 



H, = 



V 



I Br 

-dr + - + — 

r 2 



/ + 1 Br 



-TO 



The details of further calculation depend on the sign of B and below we consider different cases separately. 



4.1.1 B >0 

First we introduce two families of solutions of the Dirac equation without AB field, — E)^ = 0: 



/ 



■ -jBjl . 



■ +VB/2 ; 



Cz'e-'/^l + ^B/2 



2B \ 



(4.4) 



(4.5) 



with 



C+ = 



vn — tj\ 



2B\ 



±1/4 



(4.6) 



Define ^± (z, ( 



*±(z 



z, ( 



as in (|2.3ip . The functions ^'+(2,6') and *I'_(z,6') are delimited 



by the horizontal branch cuts I — 00 + j(95 + tt + 27rZ), In I \/ B /2 r \ +i{Lp + n + 27rZ) in the 6'-plane, whereas 



the branch cuts for *-(z, 6) and *+(z, 6^) are - In ( \jB/2 r ) + i((/3 + tt + 27rZ),cx) + 1(93 + tt + 27rZ)) (see 



Fig. 5). The arguments of both 1 + \/ B /2 ze ^ and 1 + i/_B/2 z are fixed to be zero on the line Im f ~ tt. 
Next introduce the following functions: 



w['\z) = 



Co{z) 







C±{z) 



= 27ri e 



g*('+l)v 

j*''^ 







=*(i+i)v 



wj^)(z) 



Co(z) 



+ 



■^/'•^\z)=T / e-('+i/2)«^+(z,0)d^ = 2^ze 



— inl 



-ilif 









(4.7) 




(4.8) 




(4.9) 




(4.10) 


. 5. Notice that 


(11^+) f \ 
wl '{z) 



and w|^^'^''(z) are well-defined for I > —1, while w|^^' \z), w|^^' ''(z) are well-defined for I < 0. The 



functions i/7j^''(r), ^[^^'^■'(r) satisfy radial Dirac equation (Hi — E)wi = 0. They are explicitly given by 



X 



/2B \2B 



A2 



A 



/2B \2B 



A2 



-Br^/4 



C-\^r)-'u{^ + 1,1-1,^) 



C^{^2ry'-'U[^^^1,^] 



C+{^/Bj2r)'+^U [^ + 1 + 1,2 + 1,^ 



(4.11) 
(4.12) 
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; 




c^m-- 
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-(|)-7I 
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ln\/B/2 r 

(|)+3ji- 


1) 


(|)+57i/2 

< -^^^ 
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< 
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-ln^B/2r 



(b) 

Fig. 5: Branch cuts and integration contours in the 0-plane for 
a) 1'+(2,6l) and *_(z,6i) b) ^'_(z,6l) and *+(z,6l) 



.in.. 



r(|| + i)r(i + o 



/ 



V "r(|i)r(2 + o 



(4.13) 



(r) = 



V 



r(i-0 
r(-0 



(/b72 



(4.14) 



where M{a,(3,s) and U{a,(3,s) denote Kummer's functions, see p. The solutions (|4.1ip - (|4.14p have the 
same integrabihty properties, as in the case of the disk. Namely, 

• (r) is the only solution of the radial Dirac equation, square integrable as r — > oo with respect to 
the measure rf/i^ = rdr. 

• For I > {I < —1), 'w\^^'^\r) (resp. w^^^' \r)) is the only solution, square integrable as r ^ 0. For 



I e (—1,0) both Wj^^'*''(r) are square integrable as 



• Wj^'(r) and Wj^^'^''(r) are linearly independent for / > — 1, wY'{r) and w^^' '{r) are linearly inde- 
pendent for ^ < 0. The determinant of the fundamental matrix of solutions is given by 



,{11,+} 



(I) I 



.in,-), 



det {w[^\'w\^^''^^ 



(4.15) 
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This implies 



Proposition 4.1 ([12j) Let douiHi = C^(M+). Then for I e (-cx),-l] U [0, oo) Hi is essentially self- 
adjoint. In the case I G (—1,0) it admits a one-parameter family of SAEs. 

We will consider two particular SAEs of Hjj , whose domains are composed of functions with regular at r = 
lower (or upper) component. Similarly to the above, these two cases will be referred to as corresponding to 

the SAE parameter ^ ~ '2 ^~'2^' 

Using (|4.15p . one may show that for I e (—00, —1] U [0, 00) the radial Green function Ge i{r, r') is given 

by 



Gij,/(r,r') = 



2 



{r)^[wl'\r')) 
(r) <g) (w["'*)(r') 



for < r < r' < 1, 
for < r' < r < 1. 



(4.16) 



As in Subsection [Ol the sign "+" ("-") in corresponds to / > (^ < -1). For / e (-1,0), the 

relation l|4.16p gives the Green function of the SAE with ® = (when taken with the sign "— ") and with 

e = I (for the sign "+"). 

We can now repeat word-for-word the calculation of the full one-vortex Green function G{z, z') from 
Subsection 12. 5[ using the formula (|4.16p and contour integral representations l|4.7p - (|4.10p . Setting for defi- 
niteness ^ ~ ~'2' '^^^ obtains the following representation for the Green function of the Dirac hamiltonian 
631): 



G{z,z') = 



e-Mv-v'+2^) Q{0) _,_ ^(^^ foj. ^ _ £ (_27r, -tt) 

e-i-'iv-v') G-W (^^ _^ A(z, z') for - ip' e (-tt, tt). 



with 



G(")(z,z') 



in 



d0^-iz,e)(g>'^iiz',9), 



Co{z) 



A{z,z') = Xe-'^^f-f'^ 



1 - e- 



8«7r2 



dOi jde2 '^-{z,ei)(g)i'l{z',e2) 

Co(z) Ime2 = ¥>' 



g»l-»2 _ 1 



(4.17) 

(4.18) 
(4.19) 



The integrals (|4?T8)) and (|4A9)l reduce to 



G(°)(z,z') = eT 



V 



Cii (w(z, z')) 

-^r^C2.[u{z^z)) 



\/u{z,z') 



Cl2{u{z,z'] 
-(22(^(2,2')) 



A(z,z') = 



smTTi^ 



dO 



Cii(t;(r,r',0)^ 



g— ^ rr' sinh^ ^ 



(4.20) 



(4.21) 



l^Ci2(«(r,r',0)) 
C21 (f (r, r', 0)) e«+^(^-^')C22 {v{r, r' ,6)) 



\ 



re 



\ re" + r' 

where C(''^)i u{z, z'), and w(r, r', 0) are defined as follows: 



C(«) 



2n 



IB / A2 

r — + 1 

2 V2B 



2 '^1^25^-^'^' 2 



/^2 jj ( \^ 1 -B« 
^+■^1^25'^' 2 



(4.22) 



30 



u(z, z') |z - z'p, v{r, r', 9) ^ + r'^ + 2rr' cosh 6*. 



(4.23) 



Remark. Green function for ® ~ "^^^ t)e computed in a completely analogous manner. Final answer 



differs from (|4.17p . (|4.20p - (|4.23p in only one point: one has to replace 



by —e'^ in the first 



line of the integral representation (|4.2ip for A(z, z'), just as in (|2.49p vs. (|2.47p in the disk case. 



4.1.2 B <Q 

In the previous subsection, we have obtained G{z, z') imitating the calculation made on the Poincare disk. 
Alternatively, one can simply consider the limit 

r , r' 

R oo, r —> — , r — > — 
R' R 

of the relations (|2.40p . (|2.45p - (|2.48p . using asymptotic properties of hypergeometric functions. For B < 

TT 

and ® = (regular upper component) this gives a representation of the one- vortex Green function, which 
has exactly the same form as (|4.17p . (|4.20p - (l4.2ip . but with 



C(") = 



and 



e 4 
27r 



2B 



r<-Z TT I A- 1 |_B|« \ . / TT I ^' I 1 9 |J3|M 

9T7?T'^'^5~ V ^ 91731 + -l^; 



2|Sp 



\B\u 
2 



\B\u 



2\B\ 



\B\u 
2 



^ (2^01 + 1' 2; ^"+^(2751 + 



1,1, 2 



C+ = 



m ~ E 
m + E 



1/4 



2 \ 1/4 



2i? 



(4.24) 



(4.25) 



4.1.3 5 = 

Taking further hmit B ^ in (g^Ol)- (11221), we get (still for 9 = the modification for 6 
described above) 



G(")(z,z') = 



Cii(u(z,z')) 

— C21 w z,z ) 

z' — z \ 



Vu{z,z') / \ \ 

C12 u z,z ) 

z' — z \ / 

-(22(^(2, z')) 



= — IS as 



(4.26) 



A(z,z') 



sm TTi/ 



de 



,(i+u)e+i{<p-<p') 



V 



3^+'''^^i;(r, r',6l) 



re" + r 

where the matrix ^(u) is given by 



Ci2f«(r, r',0) 
C21 («(r, r', 0)) e''+'('^-'^')C22 {vir, r' , 6)) 



C(u) ^+'^° ^^^^ ^^^^ 



C4 



m — E 
m + E 



1/4 



(4.27) 



(4.28) 



and Kq^i{s) denote modified Bessel functions. 

Using integral representations for Kq^i{s), one can write (|4.27p in a different form. Namely, for any a G 

7T TT 

such that — a| < — , \tp' — a| < — we have 



A(z,z') = e-(---')^|^ 



-Xr cosh{0 1 -\-i(a—i^)) — Xr' cosh(02+'i(ck— > 
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e0i-f)2 + 1 



C4 



'?2 + " 



Example. For < (/? < tt and < yj' < tt, we may take a = — so that 



(4.29) 



— 00 —00 



As we will see a bit later, this formula makes the computation of zero-field form factors particularly simple. 
Note that similar expressions for the one-vortex Green function on the plane have already appeared in 
different papers (see, e. g., [23ll27] ). 

4.2 Two-point tau function and Painleve V 

In order to write the tau function as a Fredholm determinant, we will consider free Dirac equation in another 
gauge. Set the potential of the uniform magnetic field to be 



(4.30) 



so that the corresponding Dirac haniiltonian 



m dx ~ idy — iBy 

-dx — idy + iBy 



-m 



commutes with the ai-momentum operator Px — —idx- The eigenspace of Px with momentum p is spanned 
by the functions g{p,y)e^^^ . Let us look at the solutions of the partial Dirac equation 



{Hj,~~E)g{p,y) = Q, 



m 



-i{dy+p-By) 



-i{dy-p + By) 



-m 



(4.31) 



As above, we assume that E is real and \E\ < m. It is convenient to choose two linearly independent 
solutions of (|4.31|) as follows: 



B>0 



B <0 



B^O 



<I>(±J(p,y) 



<i>(±)(p,y) = 



^rf|^ + i 



1/2 



1/2 



C^'D_^_^[±V2B{y-%: 



I 



V2 



1± 



CtD_^^[±^\[y^^^ 

±^C^D^^^^^(±^\(y+^^ 
1/2 \ 



±iC+ 1 T 



1/2 



(4.32) 



(4.33) 



(4.34) 



Here, Da{s) denotes the parabolic cylinder function and the constant C+ in ()4.32p . ()4.33p and (I4.34|) is 
determined by (|4.6p . (I4.25|) and ()4.28p . correspondingly. Note that <i>(+^ {p, y) ($(^' {p, y)) is square integrable 
as ?/ — *■ 00 (resp. y — > —00). These two solutions satisfy symmetry relations, analogous to (13. 39^ : 



$W(p,j/) =a,$(±)(p,y). 



(4.35) 
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The normalization in (|4.32p - (|4.34p was chosen so that in all three cases 



(4.36) 



We now adapt the reasoning of Subsection l3 .3 . 1 I to flat space. Consider a line Ly{o) = {{x,y) e | ?/ = y^"^ } 
and an arbitrary C^-valued function gy(o){x) G H^^^{£y(o)), written as Fourier integral 



9y(0) (x) 



dp .g(p,y(0))e'f^ 



(4.37) 



Decompose Fourier transform g{p,y''^^) as follows: 

g{p, y^"^ ) = 9+ {P, {p, y^"^ ) + g- {p, y^") )$(-) {p, ), 

where y'^"^) denote the functions defined by (|4.32|) . (|4.33p or (|4.34p . depending on the value of B. 

The formula (|4.36p and symmetry relations (|4.35p imply that 



9± 



{p,y^°^) = (<i>(^)(p,2/^"^))^'^. 9{p,y^"^)- 



(4.38) 



Recall that g+{p, y^^^) and g-{p, y^*^^) can be thought of as coordinates in the spaces of boundary values of 
solutions of the free Dirac equation (-fft^"* — E)ip = in the half planes y > and y < 
It is now straightforward to write down the analogs of the Propositions 13.81 and [ 



Proposition 4.2 Let us consider a strip 5 = |(a;,y) G 
can be continued to S as a solution of the free Dirac equation ( H 



AO) 



< y^^^y Suppose that V' G H^/^{dS) 



E]^|J = 0. Then 



^L.-ip) 

4'R,+ {p) 



1 

1 



'0L,+ (P) 
■4^F!.-{P) 



(4.39) 



Proposition 4.3 Assume that the strip S contains one branching point oq (i. e. y^^^ < aoy < y'^^^ ) and 
introduce a horizontal branch cut I — (— oo + iaoy, + iaoy] . Suppose that ■0 g H^^^{dS) is the boundary 
value of a multivalued solution of the free Dirac equation on S\£, which is characterized by the monodromy 



whe 



at the point oq- Then 



IpL-ip) 
V'ii,+ (P) 



(as{ao)iJjL.+^ (p) 
[psiao)i^R,-) (p) 
(75(ao)V'L,+) (p) 
('5s(ao)V'i?-) (p) 



75 (ao) Ss{ao) 



'<pL,+ ip) 
•tpR-ip) 



A^^^'^Hp^q) ^L,- 


\-{q) dq, 


(4.40) 


G^_^"'''\p,q) i,n,^ 


_(g) dq, 


(4.41) 


G^^°--\p,q) i^L,, 


-{q) dq, 


(4.42) 


A^:"'-\p,q) i^R,. 


-(<?) dq. 


(4.43) 



A, 



q) 



1 

2^ 



dxdx' ^^P-+^^-' U(^){p,y))\, At''''\z,z') a,$(T)(g,y'), (4.44) 



-oo — oo 
oo oo 



1 

'2^ 



dxdx' ($(±)(p,y)y Gt"'''\z,z') 



y'^aoy 



v,$(T)(g,y')- (4.45) 
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Here, G'(°''''^^(z, z') denotes the Green Junction of the Dirac hamiltonian h[^^ on the plane with one branching 
point ao, and A^""'"' (z, z') = G^"°^")(z, 2') - GfJiz,z'). 

The definition of tlie tau function of the Dirac hamiltonian on the plane with two branch points ai and 
a2 is also completely analogous to the disk case. Repeating the arguments of Subsection 13.11 one ends up 
with the following Fredholm determinant representation: 

r(a) = det (l - a{a2)S{ai)) , (4.46) 

where 0(02) and 6{ai) are given by (I4.40p . (I4.43|) . As above, the fact that the tau function depends only on 
the distance between the points ai and a2 allows us to choose ai = 0, a2 — t + iO {t £ and the invariance 
of with respect to x-translations reduces the problem of calculation of r(a) to finding the form factors 
A'£''^\p,q). Finally, the symmetry of the free Dirac hamiltonian iJ^^-* combined with the relations (|4.35p 
implies that 



Af-\p,q) = Af''\p,q) = Af-'\q,p) ^ A^'-\^p,-q). (4.47) 



The form factors Aj_ '"^ {p, q) are determined by the relation (|4.44p or by the equivalent formula 



1 

2^ 



dx dx' e-*P^+'"?-' A^O'"^) (z, z') = A^"'") (p, q) (p, y) ® (g, y')) ^ ■ (4-48) 



—00 —00 



y,y'^o 



We remark that Aj"'"^ (z, z') and the function A(z, z') defined by (for B ^ 0) or by (|i:7f|) - (|i:^ 

(for B = 0) are related by 

Ai°'"^(z, z') = e-(^-^') X eT^("?'-"V) ^ A(z, z'), ^ e (-^, tt). (4.49) 

Here, the first factor corresponds to a singular gauge transformation removing the AB field, and the second 
one comes from the change of the vector potential of the uniform magnetic field from (|4.ip to (|4.30p . 

Also note that y and y' in (|4.44p and ()4.48p can be chosen arbitrarily. Analogous observation in the disk 
case allowed us to obtain a more explicit representation for Aj?''^^(j>, g) by analyzing the asymptotics of a 
relation similar to (|4.48p near the disk boundary. For B the asymptotic analysis of the LHS of (|4.48p 
as y, y' —^ ±00 becomes rather complicated and we have not managed to repeat the above trick in this case. 
However, when both p, q are positive or negative, one can choose y and y' in such a way that the arguments 
of parabolic cylinder functions in the RHS of one of the relations (|4.48p are equal to zero. This leads to a 

simpler (than (I4.48P for general y,y') representation of A^2'''\pt1)- 

7T p q 

Example. Assume that B>0, p>0, q>0 and 8 = . Then, setting in (|4.48p y = — , y' = — and 

2 B B 

taking into account (|4.49p . one finds a triple integral representation for A^^''^\p, q): 

2 



2l+A''/2B 



smTTi/ / ,„ / , , I e 



Af''\p,q) = V2B-^ y^-^ — I de I d^ I V 



a: 

45 



52 (27r)5/2 TT J J ^ ^e+^{v-fi') + l 

-00 

I 1 / p2 g2 2pqe^ „ n „ n , 

2 -2 , exp<^-— — ^+ . 2 + +^^P Ctgip~2iq Ctgifi 



sin (fi sin (fi' [ 4_B \sin sin ip' sin ip sin ip' 



2B 2B \sin iy9 sin 1^9' sin 93 sin 93 

Much simpler results can be obtained for B = 0. In this case, it is convenient to introduce instead of the 
momentum p a rapidity variable 6p defined by 



A sinh 9p , \/X^ + p^ = A cosh Op . 
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Partial waves (|4.34|) can then be written as 



TT 

Set ^ ~ ~'2 substitute the representation (|4.29p for A(z,z') (recall that it is valid for y,y' > 0) into 
(|4.48|) - (|4.49p . After elementary integration over x and x' in the LHS of (|4.48p we find 

a(o.). ^ 1 sin^i. eT(i+-)(»p+«.) ^^^^^ 

AV {p, q) = = s— 3— . 4.50 

^ '"^'^^ A^coshflpcosh^g TT 2cosh^ 

TT 

Similarly, for Q ~ O'^s obtains 

sinTTzy e=F'^(^p+''9) 



A^cosh^p cosh 9q tt 2 cosh 



Ar^b^?) = . , ■ (4-51) 



It should be mentioned that the formulas equivalent to (|4.50p - (l4.5ip were first obtained by Schroer and 
Truong [32] • In a context similar to ours, they were rediscovered by Palmer in [27j . 

We finally comment on the limiting form of the equation (|l.ip in fiat space. Introduce s = R~^t and let 
R oo. Then, setting w{s) = 1 — y(t) and using asymptotic behaviour of the parameters (3 and S in (jl.3p 

^^^rr^-E^ + B^l+j^^^ ^' = , (4.52) 
it is straightforward to check that y{t) satisfies Painleve V equation 

d^y _ f I , 1 \ fdyV \dy (y-lf ( , , /3'\ , 7'y , (5'y(y + 1) 



dt^ \2y y-lJ\dtJ tdf v"'"^7j+~ + ~7^ ^'-'^^ 

with parameters a' = a = — , /?' = and 7', (5' defined by (I4.52p . The equation (|1.2p in the planar limit 
transforms into 

d f dyY X^y r^O ty t^y 

'dt'^'^^'^^^^W^A'dt) - — + T^"T(^3T)i' (^-'4) 

where rj — — — and rj{9 + 1) = 7'. The RHS of (|4.54p coincides, up to addition of a constant, with the 

Okamoto hamiltonian for Painleve V equation (|4.53p and thus the r-function (|4.46p associated to the Dirac 
operator in flat space is very simply related to Painleve V r-function. Although one could expect a similar 
relation in the case of Painleve VI and Dirac operator on the hyperbolic disk, this appears not to be the 
case l29l. 



5 Concluding remarks 

An important problem left outside the scope of this paper concerns the short-distance (s —> 0) behaviour of 
the PBT r-function. Extending the conjecture of [7] to the case 6 7^ 0, one could assume that r(s) coincides 
with the two-point correlator of twist fields in the Dirac theory in a more general classical background 
(Poincare metric -f uniform magnetic field) . External fields drastically change the infrared asymptotics of the 
correlation function, but they should not affect the exponent a in its conformal behaviour r(s — > 0) ~ Cs'' . 
To prove this rigorously, one would require a generalization of the connection formulas for Painleve VI '16[|17| 
to non-generic values of parameters (recall that in our case 7 = 0). These formulas, however, are not known 
except for the special case of PVI with /3 = 7 = 0, [fT| [T5]. 
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The theory of Painleve equations does not provide any answer for the coefficient C . For the r-function 
arising in the scahng hmit of the two-dimensional Ising model, this constant was extracted from a careful 
asymptotic analysis of the corresponding Fredholm determinant [33] . Although it seems hopeless to repeat 
such an analysis with the determinant (jl.4[) . one could try to obtain C using QFT arguments: it may be 
expressed in terms of the vacuum expectation values of twist fields, which can be computed by the method 
of angular quantization [71 [22] . 

It is curious to note that another one-parameter class of solutions of the Painleve VI equation with one 
singular parameter arises in the representation theory of the infinite-dimensional unitary group 3J. However, 
its relation to PVI transcendents studied in the present paper remains rather obscure. 



Appendix 

We wish to show that the logarithmic derivative of the r-function is invariant under the joint SU{\,1)- 
transformation of the branch point positions: 

aa.j -1-/3 f a (3 \ 
paj + a \ P Oi J 

It is sufficient to prove that the 1-form ()3.88|) is annihilated by two vector fields 

j 

whose integral curves are orbits of a compact and a noncompact subgroup of S'C/(1, 1). 

Let us first apply X to dlnr. Substituting (j3.83|) - (l3.84p into the resulting expression, one finds after 
some simplification 

(A.l) 

Two sums in (|A.ip are separately equal to zero. For example, the first one can be rewritten as 



\ Tr{[ai, _ [di, = i Tr{(e - A + hlf -{h-A + blf} = 
^TT[fg-gf-2{A-bl)ie^h)} 



= 0. 

Besides the relations (13.73^ and (|3.79p . in the above we have used the fact that the diagonal parts of the 
commutators [ai. A] and [di. A] and of the difference e — h are equal to zero. Similarly, the second sum in 
(|A.1|1 gives 

Trf/iciA + gAhiA) = Trf (AbiA - hi)AciA + (ci - AciA)AhiA) = 0. 



Next consider the action of Y. We get 
Yidlnr) = J2 El^J'^ - - f"^" - ^'))k,}+ E(«^- + ^M+i^j - 1/2, b) + 

+ E{-/'jfe(^^i - - gjkiAhi - AbiA2)fcj}. (A.2) 
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The first sum in (|A.2[) can be transformed into 

= -Tr (^e^A + h^ A- 2{A- 61) (eA + hA) + {A - blf{A + A) 
= -TT(e^A + h^A-{A-blf{A + A) 



while the third one gives 

Tr(/(lci - c,A + gA) - g{Ah, - h,A ~ /A)) - Tr(/gA + gfA). 

Summing up the three contributions in (|A.2p and using (|3.79p once again, one finds tliat rflnr is invariant 
under the flow of Y. 
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